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Abstract. The radiative cooling time of the X-ray-emitting plasma near the center in many clusters of galaxies is 
shorter than the age of the cluster, but neither the expected large drop in central temperature -nor the expected 
mass flow towards the pressure-depleted cluster centers- are observed. We propose a solution to these "cooling 
flow" problems wherein energy is supplied to the plasma by hadronic cosmic rays (CRs). The solution requires 
an energy deposition more intense and more distributed than in conventional CR models, but this alternative 
model is also successful in describing the origin of CRs, as well as the properties of gamma ray bursts and of 
the diffuse gamma background radiation. The X-ray energy emitted by clusters is supplied, in a quasi-steady 
state, by the hadronic CRs, which act as "warming rays" . The temperature distribution in the intracluster space 
is successfully predicted from the measured plasma-density distribution. Four other puzzling features of clusters 
can also be explained in simple terms: the discrepancy between their "virial" and "lensing" masses, their large 
magnetic fields, the correlation between their optical and X-ray luminosities, and the non-thermal tail of their 
X-ray spectrum. 
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1. Introduction 

The intergalactic plasma in the central regions of many 
clusters is cooler than in their outskirts, and it radiates 
X-rays at such a rate that the plasma cooling time is 
much shorter than the age of the cluster. It has been ar- 
gued (Cowie & Binney 1997; Fabian & Nulsen 1977; for 
reviews see, e.g. Sarazin 1988; Fabian 1994, 2002; Binney 
2001) that in such clusters there should be a flow of plasma 
from the outer hot regions to the cooler center — a "cool- 
ing flow" (CF) — in order to maintain hydrostatic equi- 
librium. For decades, the flows have been increasingly 
at variance with observation. Two of the old puzzles are 
that the inflowing material disappears with no observed 
trace and that the central-temperature depression is not 
as deep as expected on the basis of the plasma's cool- 
ing rate, with the observed central temperature Ti nner 
settling at a fraction (or the order of 1/2) of the outer 
temperature T outer (see, e.g. McNamara 1997). More re- 
cently, high spatial resolution imaging with Chandra (e.g. 
McNamara et al. 2000; Fabian 2000; Blanton et al. 2001; 
Allen, Schmidt & Fabian 2001; Lewis, Stocke & Buote 
2002; Blanton et al. 2002) and high spectral resolution 
measurements with the Reflecting Grating Spectrometer 
of XMM- Newton (Peterson et al. 2001, Tamura et al. 2001, 



Kaastra et al. 2001; Kahn et al. 2002; Peterson et al. 
2002) and with XMM/EPIC (Bohringer et al. 2001a,b,c; 
Molendi and Pizzolato 2001; Matsushita et al. 2002) have 
sharpened the problem: lines corresponding to gas below 
Tinner are not seen, as if the gas cooled to this temper- 
ature and then vanished. We refer to these as the "CF 
problems" . 

It has been recognized that some form of heating, e.g. 
by a central radio source, may alleviate the CF problems 
(Bohringer et al. 2002; Churazov et al. 2002). One diffi- 
culty with solutions along these lines is that the heating 
has to be distributed, like the plasma itself is. Another 
difficulty is that the bremsstrahlung cooling rate is pro- 
portional to r? e T 1 / 2 , with n e the plasma's electron number 
density: a suitable heating and pressure-building mecha- 
nism must somehow adapt itself to this behaviour, par- 
ticularly as a function of n e , which varies by orders of 
magnitude along the clusters' radii 1 . 

Heating by cosmic rays (CRs) has also been discussed 
(Rephaeli 1987; Rephaeli & Silk 1995), but it does not 
seem to be a good solution. The traditional theory of CRs 
posits that the bulk of them is accelerated by expanding 
supernova (SN) shells. The total energy in CRs made this 
way falls short of the required heating energy by more 
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1 Hydrogen in the plasma is fully ionized, and we are using 
in the Introduction n p ~ n e . 
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than one order of magnitude. Moreover, the acceleration 
takes place close to the central parts of galaxies, and the 
CRs would not deposit energy throughout a cluster as 
uniformly as it is required. 

A somewhat different theory of the origin of CRs has 
been proposed, according to which SNe indeed accelerate 
CRs, but in an unconventional way. Like matter-accreting 
quasars and microquasars, most SN explosions would be 
accompanied by the emission of bipolar jets of ultrarela- 
tivistic "cannonballs" (CBs), with initial Lorentz factors, 
70, of 0(1O 3 ). The asymmetry in these jets would be re- 
sponsible for the observed large "kick" velocity of neu- 
tron stars (Dar & Plaga 1999; Dar & De Rujula 2000a). 
When these CBs are moving close to the line of sight of 
an observer, they are seen as gamma ray bursts (GRBs) 
and their afterglows, providing an extremely simple, de- 
tailed and successful description of these phenomena (Dar 
& De Rujula 2000b; Dado, Dar & De Rujula 2002, 2003a, 
2003b; De Rujula 2002; Dar 2003 and references therein). 
The CBs travel for kiloparsecs in the interstellar medium, 
decelerating by knocking out its constituents, which are 
thereby accelerated to become CRs. This CR theory pro- 
vides a qualitative understanding of the CR spectrum 
(Dar & Plaga 1999; Plaga 2002) and a simple and very 
satisfactory understanding of the diffuse, high-galactic- 
altitude Gamma Background Radiation (GBR), which 
originates in the halo of our Galaxy, and is therefore not 
cosmological (Dar & De Rujula 2001a). The CRs in this 
theory are not confined to the disk of the Galaxy, but per- 
meate a much larger halo. The theory also implies a total 
CR production rate that is more than an order of mag- 
nitude larger than the conventional figure, and yet, not 
in disagreement with observations (Dar and De Rujula 
2001b). 

We show here that this unconventional theory of CR 
production resolves in a very simple way the CF prob- 
lems and results in a predicted temperature distribution 
in clusters that is in very good agreement with the obser- 
vations, while implying the absence of significant flows. 
What are needed, we posit, are not cooling flows, but 
"warming rays" (WRs): cosmic rays 2 that supply the X- 
ray energy emitted by clusters from their central regions 
and play a role in determining the temperature and pres- 
sure distribution in the whole intracluster space. 

As a bonus, our theory offers plausible ways out from 
four other cluster conundra: the X-ray/optical luminos- 
ity correlation, the lcnsing/X-ray mass discrepancy, the 
(quite) large intracluster magnetic fields, and the origin 
of the non-thermal tail in the hard X-ray spectrum: 

— The luminosities of clusters in optical and X-ray bands 
are observed to be correlated: they are proportional 
(e.g. Miller & Nichol 2001). Once again, this requires 
an apparently surprising concatenation between two 

2 We do not use the expression "warming cosmic rays" , since 
cosmic rays are not "cosmic" in the current sense of the word. 
When referring to CRs in CF clusters, we use the expressions 
"CRs" and "WRs" as synonyms. 



very different emission mechanisms. This correlation, 
and the value of its proportionality factor, are the ones 
implied by our theory of WRs. 

— The mass in cluster cores estimated from X-ray mea- 
surements is systematically a few times smaller than 
that estimated from observations of gravitational lens- 
ing (e.g. Hattori et al. 1999, and references therein). 
The X-ray estimate relies on an approximate hydro- 
static equilibrium of the plasma in the clusters. If 
magnetic fields, WRs and the plasma they heat up 
are in energy equipartition, the pressure gradient of 
the first two not taken into account in the assumed 
equilibrium — is of the right order of magnitude to 
solve the mass-discrepancy problem. 

— The very large CR production rate that we pre- 
suppose corresponds — again for an assumed energy 
equipartition — to the surprisingly large fields that are 
observed in clusters (e.g. Carilli & Taylor 2002; Eilek 
& Owen 2002, and references therein). 

— The high energy tail of the spectrum of the X-ray emis- 
sion observed in some nearby clusters is non-thermal: 
dn 7 /dE ~ E- 2 (Fusco-Fcmianno et al. 1999, 2000; 
Rephaeli, Gruber & Blanco 1999). This is the expec- 
tation for the emission from the knocked on electrons 
produced in the CR showers that, in our theory, "warm 
up" the inner-cluster plasma. 

We summarize in the text the features of the CB model 
most relevant to the problem of CF clusters, relegating to 
an Appendix the outline of their derivation, as well as a 
short overview of the model itself. 

2. The distribution and luminosity of CRs 

We have argued elsewhere that the "natal kicks" of neu- 
tron stars, cosmic rays, long-duration gamma-ray bursts 
and the diffuse gamma background radiation have a com- 
mon origin: highly relativistic jets of cannonballs emit- 
ted in supernova explosions (for CRs and neutron stars, 
see Dar & Plaga 1999; for the GRBs and their associ- 
ated SNe see, e.g., Dado et al. 2002; De Rujula 2002; Dar 
2003 and references therein; for the GBR, see Dar & De 
Rujula 2001a). In our unconventional view of every one 
of these subjects, high-energy astrophysical phenomena 
would have a remarkable unity. 

Of the quoted items, the most relevant to the current 
work are those leading to the estimate of the relation be- 
tween the optical luminosity and the CR luminosity of a 
galaxy, and to the distribution of CRs in space. We sum- 
marize them here. 

In Dar & Plaga 1999, Dar & De Rujula 2001a,b, and 
in the Appendix we argue that the CR luminosity of the 
Galaxy should be: 

i CH ~5x 10 42 erg s"\ (1) 

which is more than one order of magnitude bigger than 
conventionally believed. We also contend that the exten- 
sion of the CR distribution in the Galaxy may be much 
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larger than generally accepted, so that the CRs permeate 
a "cosmic-ray halo" extending for tens of kiloparsecs. Such 
a CR distribution in the galaxies of a dense cluster would 
permeate the cluster fairly uniformly and result in a CR 
population proportional to the number density of galax- 
ies in the cluster 3 . For the O(10) fiG magnetic fields of 
an intracluster medium (e.g. Carilli & Taylor 2002), the 
bulk of the CRs does not escape the cluster over times 
comparable to its age (e.g. Colafrancesco & Blasi 1998). 

Our arguments are based on observations of long- 
duration GRBs and core-collapse SNe (of type II, lb and 
Ic). In using them to discuss the CR luminosity of clus- 
ters we ought to ponder the relative frequencies of SNe 
of various types in spirals or ellipticals, and on the abun- 
dances of the various galaxy types in clusters. Though 
such considerations would not affect our results by large 
factors, we make the simplifying assumption that all SN 
types produce jets and CRs of similar energy: Type la su- 
pernova explosions occur as the mass of a white dwarf in a 
binary increases beyond the Chandrasekhar limit because 
of the mass accretion from a non-degenerate companion 
or a merger with another white dwarf. These processes 
result in an unstable accretion configuration akin to that 
observed in micro-quasars, and should also result in the 
emission of jets. If, as we presume, Type la SNe are asso- 
ciated with short- duration GRBs, the energy of their jets 
should not be much smaller than for core-collapse SNe. 

The contention that long-duration GRBs are produced 
by core-collapse SNe (Dar 1999, Dar & Plaga 1999, Dar 
and De Rujula 2000, Dado et al. 2002) has recently re- 
ceived striking observational support: the predictions of 
Dado, Dar& De Rujula 2003c of a SN associated with 
GRB 030329, its properties and even the date when it 
would be convincingly discovered have been found to be 
correct (Stanek et al. 2003). 



3. Cosmic rays in the intracluster (IC) medium 

The IC space of dense clusters is filled with a highly ion- 
ized plasma at a temperature T of a few keV, electron 
number density n e = 10 -2 to 10 -3 cm -3 near the center, 
and slightly sub-solar metal abundance. This plasma is re- 
sponsible for most of the cluster's thermal X-ray emission, 
via bremsstrahlung in the 0.1 to 10 keV range (see, e.g. 
Sarazin 1988). 

In the CB model the main acceleration mechanism of 
electrons and nuclei by CBs is simply the "collisionless" 
scattering of these ambient particles by the magnetic field 
of the moving CBs, which makes them recoil with energies 
proportional to their masses. Thus, the ratio of electron 
to nuclear CR energies is ~ m e /M, with M the mass of 
the corresponding nucleus: the bulk of the CR energy is 
carried by CR nuclei, mainly protons. 



3 Active Galactic Nuclei (AGNs) emit powerful jets of CBs 
that should also be efficient CR accelerators, but not all clus- 
ters contain observed AGNs, nor is the CR population that 
AGNs generate uniformly distributed. 



The dominant energy-loss mechanisms for hadronic 
CRs in the IC medium are Coulomb collisions with the am- 
bient electrons, which lead to electromagnetic showers ini- 
tiated by the knocked-on electrons, and hadronic collisions 
with the ambient nuclei, which also lead to high-energy 
electrons and positrons via pion production, followed by 
the decay chains ir + — > /i + v, tt~ — > /i~ D; ^ — > e^vv. 
A relatively small fraction of the available energy escapes, 
carried away by neutrinos and by 7 rays from ir° produc- 
tion and decay. The rest of the CR energy is deposited in 
the plasma and thcrmalizes, via electromagnetic interac- 
tions of the electrons and positrons, and elastic collisions 
of the low energy hadronic CRs with the IC electrons and 
nuclei. 

A highly ionized plasma is much more efficient than 
non-ionized gas in slowing down the CR nuclei. Energy 
losses via Coulomb collisions with neutral atoms become 
significant only for energy transfers larger than their ion- 
ization potential I. The corresponding small impact pa- 
rameters imply relatively small cross sections. For the IC 
plasma the role of I is played by the much smaller quan- 
tity I p = hujp, with ujp = [47r n e e 2 /m e } 1 ^ 2 the plasma fre- 
quency. Let Z 2 be the (suitably averaged) squared charge 
of the plasma's nuclei. For moderately relativistic CR nu- 
clei with a charge z moving with a velocity v = (5 c and 
Lorentz factor T, the energy-loss rate is: 

: 2 [3 2 T 2 



dE 
~dt 



Kz z Z l 



1 





In- 



2 m P 



with K = 4 7r rip r 2 



m e c 



and 



2 /m e 



(2) 



2.82 fm. 



The rate of Eq. (J2J , for typical plasma densities near the 
center of clusters (n e ~ 10~ 2 cm -3 ), is ~ 4 times larger 
than for neutral hydrogen. For protons (z = 1) with a 
typical CR energy (T — 2) and for a plasma consisting 
of - 0.93% H, - 0.07% He (by number), and traces of 
heavier elements, Z 2 ~ 1.14 and dE/dt » 10~ 17 GeV/s, 
for the quoted electron density. 

Hadronic collisions between CRs and plasma nuclei 
constitute an energy-loss mechanism approximately as ef- 
ficient as that of Eq. Consider the collisions of CR 
protons with ambient protons, the discussion of proton- 
nucleon and nucleon-nucleon collisions being far too cum- 
bersome to be justified by the small admixtures of nuclei 
heavier than protons in the CR and IC constituencies. For 
r > 2 of the incident proton the pp cross section is o~ pp ~ 40 
mbarn, and dominantly inelastic. In practically all of these 
collisions the incident proton survives unscathed, as the 
"leading" (most energetic) final-state particle, carrying on 
average some 70% of the incident proton's energy. Thus, 
a proton's energy-loss rate by hadronic collisions in a hy- 
drogenic plasma is: 



— ~ 0.7 n K a. 
at 



i>i> 



cE: 



(3) 



that is dE/dt ~ 1.6 x 10~ 17 GeV/s, for our reference n e 
and r = 2. This is very close to the Coulomb energy loss 
computed at the end of the previous paragraph for the 
same reference parameters. 
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The characteristic time for a CR proton to lose its 
energy is r = Ej (dE/dt), with dE/dt the sum of Eqs. (01 
and Once again for the chosen reference values, r ~ 
4 x 10 16 s. This is shorter than the age of the cluster 
which, for the relatively nearby clusters of interest here, 
is approximately the age of the Universe, fy~3x 10 17 s. 

To summarize, to a good approximation, CRs deposit 
their energy within a CF cluster in a time shorter than 
the cluster's age. Thus, in the absence of other significant 
energy sources, the IC plasma cools until its energy-loss 
rate by bremsstrahlung is balanced by the CR energy- 
deposition rate. Since the source of X-ray energy at the 
cluster centers are the "warming" CRs, the cluster lu- 
minosities in CRs and X-rays are approximately equal: 



4. The X-ray luminosity of a cluster 

"Cooling Flow" clusters contain large amounts of plasma 
and, as we have argued, we may estimate their X-ray lumi- 
nosity as the CR luminosity of the galaxies they contain. 
The latter, in turn, can be estimated from their optical 
luminosity, L*, as follows. 

The Shapley-Ames "fiducial" sample of 342 galaxies 
within the Virgo circle has a mean-galaxy B-band lumi- 
nosity L, G) — 6.7 x 10 9 £ Q (-B) and, for a galaxy or clus- 
ter of B-band luminosity L, a mean rate of 2.07 L/ L {G) 
SNe per century (Van den Bergh and Tamman 1991). We 
have argued that the total CR energy produced by a SN 
is roughly twice the jet energy estimated from the typ- 
ical neutron-star velocities; see the Appendix, Eq. (TTSJ). 
Assembling these facts and converting from a B-band to 
an optical solar reference, we conclude that the X-ray, CR 
and optical luminosities of clusters ought to satisfy: 



0.02 L*. 



(4) 



The predicted relation between L x and is consistent 
with the observations for CF clusters, as we show in 
Fig. JTJ), adapted from Miller et al. 1999. The panel (c) 
in this figure disagrees with the prediction, but it refers 
to groups and not clusters of galaxies. Since WRs are not 
confined in the magnetic field of a group, their X-ray and 
optical luminosities should be linearly related, albeit with 
a smaller coefficient than that of Eq. as observed. 
Encouraged by these results on the optical/X-ray lumi- 
nosity relation, we proceed to study its radial distribution 
in CF clusters. 

5. The temperature profile in cooling flow clusters 

The "virial" expectation for the temperature of the plasma 
in a cluster, for a spherical distribution with mass M(r) 
within a radius r, is: 



kT(r) 



Gm p M{r) 



3r 



2keV 



M(r) 



10 13 M, 







100 kpc 



(5) 



The observed "outer" temperature of the plasma in a 
cluster, T out er, is — at large radii where the cooling time 
should be longer than the cluster's age — roughly compat- 
ible with the virial expectation. In the inner regions the 
temperature is somewhat smaller (~ 1/2 of T outer ) than 
that of Eq. JSJ, and significantly larger than the cooling 
rate would imply (see, e.g., Peterson et al. 2001). 

What is the effect of the "warming rays" on a clus- 
ter's temperature distribution? To answer this question 
we make the following simplifying assumptions: 

— Spherical symmetry. 

— The density of the plasma is proportional to that of 
galaxies and of the SNe they contain (clearly without 
sufficient "resolution" to individuate single galaxies). 
This implies that the CR distribution — which in our 
model traces with bad resolution the SN distribution — 
also traces the plasma distribution: the radial distribu- 
tion of the CR intensity is proportional to that of the 
plasma density. 

— The temperature of the plasma in its "initial state" (as 
the cluster is born) can be approximated by a radius- 
independent constant. 

— Cooling by X-ray emission and heating and pressure- 
building by WRs are the dominant evolutionary 
agents. Flows and cluster ageing are relatively unim- 
portant. 

Given our assumptions, the rate per unit volume at 
which the plasma loses energy by the emission of X-rays 
is: 

de 
dt 



= a [n e (r)] 2 y/T(r,t); 



2 n 7r 3 e\ 



4.8 x 10~ 24 



3 3 h ml c 3 

1 erg cm 3 



(6) 



where z is an average charge of the IC matter and in the 
numerics we have approximated the Gaunt factor g by 
unity. The rate per unit volume at which the WRs deposit 
energy in the plasma is also quadratic in n e : 



de 
dt 



b[n e {r)}\ 



(7) 



since we have first argued and then assumed that the WR 
production rate per unit volume is proportional to n e . The 
evolution of the temperature anywhere in the IC plasma 
satisfies the energy-conservation relation: 



3 k n e (r) 



dT(r, t) de 



dt 



dt 



de 
dt 



(8) 



with the initial condition T(r, 0) = Tj. Rewrite Eq. JSJl as: 



dT(r,t) an e (r) 
3fc 



T~b/a 



dt 



(9) 
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and integrate, to obtain: 



In 



y/T(r,t)/T f -l 




6n e (r)kTf 



(10) 



where Tf = (b/a) 2 is the steady temperature first reached 
at the center of the cluster, where the density is high, and 
at asymptotically late time in the periphery, where the 
density is vanishingly small. We have written the r.h.s. of 
Eq. Ijl U| l in a way that shows that the evolution of T(r, t) 
is determined by the ratio between the WR deposition 
rate bn 2 and the thermal energy per unit volume at the 
cluster's central temperature. 

The temperature evolution described by Eq. (fTTjl) is 
shown in Fig. @, for t — 0,t v /3,t u and 3t u , with t v the 
current age of the Universe 4 . In the particular example in 
Fig. J2J the central density is high enough for the equi- 
librium between X-ray cooling and WR heating to have 
been reached at the center of the cluster in a time of the 
order of its current age. In such a case, Tf ~ Ti nner , the 
presently observed inner-cluster temperature. We shall see 
that some clusters have reached this stage, while others 
have not. For all the clusters we study, the cooling time 
in the outer parts of the cluster is much longer than t u , 

SO that Ti ~ T outer- 



6. Temperature distributions of specific clusters 

We fit Eq. (fTTTf) to the temperature distribution of a given 
cluster as follows: we fix from the observations an ansatz 
value of Ti ~ T ou ter, we input the observed n e (r) (either 
point by point, or via a "/3-model" or a variation thereof) 
and we vary a single parameter b — or, equivalently, Tf — 
to obtain a best fit T(r, t) at the current age of the clus- 
ters (t = t(z) is fixed by the cosmological model we have 
adopted). 

We show in Figs. to JSJ the measured density pro- 
files and the predicted and observed temperature profiles 
of the clusters A1795, Sersic 159-03, Hydra A, and A2029. 
The temperature and density data are taken, respectively, 
from Tamura et al. (2001), Kaastra et al. (2001), David et 
al. (2001) and Lewis et al. (2002). These figures demon- 
strate that the correlation between the input density 
profile and the output temperature profile predicted by 
Eq. (|10f> is in excellent agreement with the observations: 
in all cases the shape of T(r) is snugly reproduced. 

The examples of A2052 and EMSS1358+6425 are 
shown in Figs. Q and (JHJ). The data are from Blanton 
et al. (2002) and Arabadjis et al. (2001), respectively. In 
the first of these clusters, the observed n e (r) is not very 
smooth and we have made two predictions for T(r). One 



4 We use a cosmology with Qa = 0.7, £Im = 0.3 and Ho = 
75 km s _1 Mpc -1 but in relating the radial coordinate of a 
given cluster to the observed angular coordinate, we use in 
each case the value of Ho assumed by the observers. 



prediction uses, as in all other cases, the best fit n e (r) 
and results in a smooth prediction for T(r). The other 
prediction uses the measured point-by-point density, and 
results in a point-by-point prediction for T(r), shown in 
Fig.Q as a series of triangles. The latter prediction follows 
the shape of the observed T(r) somewhat better than the 
smooth curve does. EMSS1358+6425 is special, in that its 
measured T(r) is compatible with a constant, albeit with 
very large errors. 

The values of the central asymptotic temperature, Tf, 
resulting from our fits are shown in Table I, and compared 
with Ti nner , the innermost observed temperature. Also 
given in the table are the values of Ti, the initial constant 
input temperature, and of T outer , the outermost observed 
temperature 5 . The meaning of all these temperatures is 
illustrated in the lower panel of Fig. J2Jl. In Sersic 159- 
03 and in Hydra, Tf — Ti nner within errors, so that the 
centers of the clusters have reached a nearly steady final 
temperature. For A1795 and A2052 there is an indication 
that Tf <Tinner, so that the cluster's central temperature 
should still cool in the future (not a prediction the success 
of which we are likely to witness). In these four clusters, 
Touter < Ti , indicating that their outer regions have also 
cooled with time. The fits to the temperature profile in 
Figs. Q to © are very satisfactory. In the case of A2052 
the observed density is so erratic that the formal errors 
of our simple fit cannot be trusted; the large % 2 reflects 
this, as is clear in Fig. J7J). For EMSS1358+6425, shown 
in Fig. |JS}, the errors in the observed temperatures are so 
large that the fit is not very illuminating. 



7. The mass discrepancy in clusters 

The gravitational lensing of distant galaxies by an inter- 
vening cluster can be used to infer the cluster's mass dis- 
tribution (e.g. Grossman & Narayan 1989). Such mass es- 
timates are quite consistent with the ones derived from 
the observed velocities of the galaxies in the clusters. On 
the other hand, the mass deduced from lensing or from the 
motion of the galaxies is ~ 3 times larger than the mass 
deduced from the clusters' X-ray emission (e.g. Hattori et 
al. 1999). Fabian and Allen (2003) obtain an agreement 
between the two mass estimates by using the density pro- 
file obtained by numerical simulations (Navarro, Frenk & 
White 1997), but this profile deviates significantly from 
the observed plasma distribution near the center (where 
it diverges as 1/r) and in the periphery (where its 1/r 3 
asymptotic behaviour is significantly steeper than the ob- 
served one). 

Assume the WRs to be in approximate energy equipar- 
tition with the magnetic field that confines them. Given 
our other assumptions on WRs, their pressure is propor- 
tional to the thermal pressure of the plasma. Let $(r) be 
the gravitational potential of a spherically symmetric clus- 



For two clusters, Sersic 159-03 and A2052, the quoted value 
of T ou t e r is that of the penultimate radial point, since the lower 
value of the last point may indicate the "end" of the cluster. 
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ter and p, the mean atomic weight of its plasma. The total 
mass-density distribution of the cluster, p tot , can be in- 
ferred from the density distribution of the X-ray-emitting 
gas, Pgas, via Poisson's equation 6 : 



Ptot 



4ttG 



V 2 $ 



(1 + a)fcr V 2 /np- 1 



(11) 



where a is the proportionality constant relating the sum 
of the magnetic and CR pressures to the thermal pressure. 
In approximate equipartition (a ~ 2), the total mass de- 
duced from the temperature and density distribution of 
the X-ray-emitting plasma is approximately three times 
larger than that inferred by neglecting the magnetic and 
CR pressures. This brings the different mass estimates to 
better agreement, as can be seen from Fig. (0, adapted 
from Hattori et al. (1999). The substantial contribution of 
the magnetic and WR pressures helps in explaining why 
our last simplifying assumption in Section 5 — that flows 
are relatively unimportant in determining the temperature 
profiles — may be a fair approximation. 

If the magnetic field is in rough energy equipartition 
with the WRs, the corresponding pressures are related by 
B 2 /(87r) — 2n e kT. The magnetic field would thus be: 



B(r) ~ (20^Gauss) 



\T(r)- 


1/2 


n e (r) 


1/2 


3keV 




lO" 2 cm- 3 





(12) 



which is in the observed domain, (e.g. Carilli and Taylor 
2002; Eilek & Owen 2002; and references therein). 

8. The X-ray tail 

CR nuclei of charge z, velocity (3 c and Lorentz factor T 
deposit their energy in the IC plasma via Coulomb scat- 
tering of the plasma electrons. These knocked-on electrons 
recoil with kinetic energies up to K max ss 2 m e c 2 T 2 . In a 
plasma of electron density n e the number of electrons per 
unit time and volume scattered to a kinetic energy K in 
the range hoj p <C K < K max is: 



d 2 N e 
dKdt 

F(K) 



2„ 2 F(K) 



2 7r r e n e m R c z 
1-/3 2 



(3K 2 ' 



(13) 



We are interested in electrons of energy much smaller than 
the K max values corresponding to typical CRs, so that the 
integration over the CR spectrum, which peaks at T— 1 = 
0(1), preserves the K~ 2 shape of Eq. i|13|) up to K = 
0(m e c 2 ). 

The scattered electrons cool and thermalize in the 
plasma mainly by bremsstrahlung emission at high en- 
ergy and by Coulomb collisions at low energy. A detailed 



6 We neglect the (logarithmic) temperature gradient relative 
to the corresponding density gradient, a good approximation 
in the cases at hand. 



calculation of the resulting quasi-steady electron energy 
distribution is beyond the scope of this paper, but this 
distribution must be thermal at low energies with a power- 
law tail whose power-law index is n ~ 2. Hence, the 
thin bremsstrahlung emission, which dominates the ra- 
diative cooling of these electrons, has a thermal (thin) 
bremsstrahlung shape at low energies, with a power-law 
tail with the same index as the energetic electrons. This 
behaviour (prior to photoabsorption corrections along the 
line of sight) may be roughly interpolated by: 



(14) 



dE 

Such "thermal" distributions with the expected power- 
law tail have been observed in the X-ray emission of 
some clusters, e.g. by BeppoSAX (e.g. Fusco-Femiano 
et al. 1999; 2000) and by RXTE (Rephaeli, Gruber & 
Blanco, 1999). Non-thermal power-law tails in the "ther- 
mal" bremsstrahlung emission from plasmas that contain 
high-energy CRs have also been observed in SN remnants 
(e.g. Mineo et al. 2001; Dyer et al. 2001), elliptical galax- 
ies (e.g., Guainnazzi & Molendi, 1999), and GRBs (e.g. 
Dar & De Rujula 2000b; De Rujula 2002; Dar 2003 and 
references therein). 

The example of the spectrum of the SN remnant SNR 
MSH 15-12 (Mineo et al. 2001) is given in Fig. lO, show- 
ing that the data can be described by Eq. (|14(l with n = 2, 
as expected, and (in keV units) T = 5 and 8 = 1.5. The fig- 
ure also shows the spectrum of the cluster A2256 (Fusco- 
Femiano et al. 2000), for which we used n — 2, T — 5 and 
5 = 3. A relatively large 5 in a given cluster is indicative of 
a large WR supply, so that its central temperature trough 
should be relatively unpronounced. There are currently no 
data to test this expected correlation. 

We conclude that the observed non-thermal power-law 
tail of the X-ray emission from CF clusters is but one 
example of the behaviour of plasmas subject to a high- 
energy CR flux, and that its presence adds consistency to 
our contention that the cluster's plasma is subject to an 
intense flux of "warming rays" . 

9. The radio emission 

The diffuse radio emission from the halo or IC medium of 
clusters was first detected by Large et al. (1959) in Coma, 
and first studied by Wilson (1970). It is due to synchrotron 
radiation in the magnetic field of the IC medium, emitted 
by electrons with Lorentz factors T e ~ 10 4±1 . A detailed 
discussion of this radio emission is beyond the scope of 
this paper, but some comments are appropriate. 

The radio-emitting electrons lose energy by inverse 
Compton scattering, with a cross-section cr T = 0.665 x 
10" 24 cm~ 2 , on the real photons of the cosmic background 
radiation and on the virtual ones of the cluster's magnetic 
field, whose energy densities are e CMB ~0.26 eV cm" 3 and 
e B = B 2 /(8ir). The electron's cooling time: 

3 m e c 



4 o T (e c 



,)r e 
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~ (7.4 x 10 15 s) 6cMB -5- , (15) 

is very much shorter than the age of the cluster. 
Consequently, these electrons must be supplied continu- 
ously and fairly uniformly over the cluster. In particular, 
they cannot be supplied by electron-accelerating SN shells, 
since the electrons' cooling time is also much shorter than 
their diffusion time out of the dense regions of galaxies. 

In our WR model, the high-energy electrons are sup- 
plied — steadily and in a distributed fashion — by the 
highly rclativistic and far-reaching CBs jetted by SNe (and 
AGNs) . The energy of these electrons is locally re-emitted 
by synchrotron radiation and dominates at radio frequen- 
cies. 

10. The 7-ray emission 

In a cluster, as we have shown, the CR protons lose about 
as much energy to knocked-on electrons as they do in 
hadronic collisions with the IC nuclei. Very roughly, the 
CR proton spectrum carries comparable kinetic energy 
above and below 1 GeV. The higher-energy protons, when 
colliding hadronically, transfer about 1/3 of their energy 
to 7r°'s, which decay into 7 rays. This implies that the 7- 
ray luminosity of a CF cluster ought to be (very) roughly 
1/6 of its X-ray luminosity. 

The 7-ray spectrum is hard to ascertain with confi- 
dence, one reason for this being that the spectral shape 
of the parent CRs in a cluster is not known. The spec- 
tral shape of the CRs on Earth is affected at low energies 
by the planet's magnetic field and, at all but the highest 
energies, by the energy-dependence of the confining time 
of CRs in the Galaxy and its halo. Predicting the spec- 
tral shape of a cluster's 7-ray spectrum, therefore, requires 
extra assumptions beyond the ones made in this paper. 

11. Conclusions 

Cooling flow clusters were named after an inference — the 
flow itself — and not a direct observation. It is therefore 
not very surprising that further observations led to contra- 
dictions. We have seen how an extremely simple model of 
the origin, luminosity and location of cosmic rays provides 
a good understanding of the main "cooling flow problem" : 
the observed temperature profile of these clusters. This 
Warming Ray model also explains in simple terms the 
clusters' X-ray luminosity, its correlation to the optical 
luminosity, the cluster's magnetic field intensity, the dis- 
crepancy between virial and lensing masses, and the high- 
energy tail of the X-ray spectra. 

The Warming Ray model has no specific consequences 
regarding other interesting cluster properties — such as 
their "metal" abundances — but it is part of a non- 
standard unified view of high-energy astrophysical phe- 
nomena, which includes the "natal kicks" of neutron stars, 
quasar and microquasar emissions, 7-ray bursts and their 
afterglows, and the origin and acceleration of cosmic rays. 



12. Appendix: Astrophysical Cannonballs 

We briefly review the CB model of GRBs, as well as our 
unconventional views on the CR distribution and luminos- 
ity of galaxies (for other reviews and for references, see e.g. 
De Riijula 2002, Plaga 2002, Dar 2003). 

12.1. The cannon 

The "cannon", or engine, producing the jets of CBs, is 
not understood. We assume that, in a SN explosion, a 
small fraction of the infalling matter may be jetted along 
the axis with highly rclativistic velocities in a succession 
of shots, before its reservoir is exhausted. This emission 
of CBs would be akin to that observed in quasars and 
microquasars, in which episodes of accretion from a disk 
or torus into a central compact object result in the axial 
ejection of relativistic "plasmoids" — or CBs — made of 
ordinary matter 7 . 

12.2. The GRB 

When crossing the SN shell and the matter distribution 
produced by the "wind" emitted by the progenitor star, 
the front of a CB is heated to a time-dependent temper- 
ature T, of O(l) keV. The typical initial Lorentz factor 
of a CB, 70 ~ 10 3 , follows from the observed energies of 
the GRB photons which, Doppler-boosted and cosmolog- 
ically red-shifted, peak at a fraction of 1 MeV. The quasi- 
thermal radiation a CB emits, boosted and collimated by 
its relativistic motion, is a single 7-ray pulse in a GRB. 
The cadence of pulses reflects the chaotic accretion and 
is not predictable, but the individual-pulse temporal and 
spectral properties are. 

12.3. The afterglow of a GRB 

The ejected CBs, as observed in micro-quasars, are as- 
sumed to contain a tangled magnetic field. As they plough 
through matter, they gather and magnetically scatter 
its baryonic constituents, mainly protons. The re-emitted 
protons exert an inward pressure on the CBs, which coun- 
ters their expansion and makes them reach an asymptotic 
radius R CB , in minutes of the observer's time 8 . The GRB's 
afterglow is dominated by synchrotron radiation by the 
electrons swept in by the CB in its voyage through the 
interstellar medium (ISM). 

12.4. The generation of CRs 

A CB ploughing through the ISM decelerates, and its 
time-dependent Lorentz factor j(t) is explicitly calculable. 

7 In the CBs of the micro-quasar SS 433, Balmer H and 
He lines (e.g. Eikenberry et al. 2001) and the Kq line of Fe 
(Migliari et al. 2002) have been observed. 

8 This mechanism may explain the surprisingly small size of 
the CBs emitted, for instance, by the Pictor-A quasar (Wilson, 
Young, & Shopbell 2001) 
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The ISM nuclei of mass M that are simply reflected by the 
CB's magnetic field at time t recoil with a flat energy dis- 
tribution extending to 2j(t) 2 Mc 2 . Their time-integrated 
spectrum is the CR spectrum, extending up to the spectral 
feature known as the "knee" . The higher-energy spectrum 
would be due to the ISM nuclei that have been Fermi- 
accelerated in the CB's tangled magnetic field before they 
exit again into the ISM. The overall CR spectrum and 
the interpretation of its very high-energy "ankle" are dis- 
cussed in Dar & Plaga 1999. The fraction of the original 
CB's energy ending up in electrons is m e /m p . Thus, 
the overall efficiency of this mechanism in accelerating CR 
nuclei is ~ 100%: essentially all the initial kinetic energy 
of a CB ends up in CR nuclei. 

12.5. The CR luminosity of the Galaxy 

What is the total energy of the CRs produced by the jets 
of CBs emitted in a single SN explosion? In the case of 
a core-collapse SN leading to the formation of a NS, we 
may roughly estimate this quantity as follows. The mean 
sky velocity of neutron stars (Lyne and Lorimer, 1964) is 
(v NS ) ~ 450 ± 90 km s _1 . Assume this velocity to be due 
to a "natal kick" : a momentum imbalance between the 
two jets of CBs emitted axially in opposite directions, as 
the NS is born. This implies: 



M 

, c ~4x 10 51 NS erg, 
1.4 Mq 



(16) 



with the inequality becoming an equality in the case of 
a large momentum asymmetry. The total energy per SN 
ending up in CRs would be of the order of twice the r.h.s. 
of Eq. Unj. l;2cl;2c The estimated rate of Type II, lb and 
Ic supernovae in the Galaxy is R SN ~ 1/50 per year (van 
den Bergh and Tammann 1991). The CR luminosity of a 
Milky- Way-equivalent galaxy would then be given by: 



2E jet R SN m 5.1 x 10 42 ergs 



(17) 



the numerical result being the one quoted in Eq. Q . This 
is more than one order of magnitude larger than the classic 
estimate L CR ~ 1.5 x 10 41 erg s _1 , based on the observed 
ratios of primary CRs to the secondary ones produced by 
collisions with the ISM (Drury et al. 1989). The alterity 
can be easily understood: the ISM gas density, volume and 
grammage used to derive the classic result all refer to a 
region close to the visible Galaxy. The derivation is invalid 
(e.g. Longair 1992) if the locally observed CRs may have 
spent a fraction of their travel-time in a much less dense, 
though magnetized halo. 

12.6. The GBR and the locus of CRs 

In Dar & De Rujula 2001a we analysed the possibility that 
the diffuse Gamma Background Radiation at high galactic 
latitudes could be dominated by inverse Compton scatter- 
ing of CR electrons on the cosmic microwave background 
radiation, and on starlight from our own galaxy. Assuming 



the mechanisms that accelerate galactic CR hadrons and 
electrons to be the same — as they are if accelerated by 
CBs — we derived simple and successful relations between 
the spectral indices of the GBR above a few MeV, and 
of the CR electrons and CR nuclei above a few GeV. 
We reproduced the observed intensity and angular depen- 
dence of the GBR, in directions away from the galactic 
disk and center, without recourse to hypothetical extra- 
galactic sources: the GBR is, we argued, predominantly 
not "cosmological" . 

The above results require that the CR population of 
the Galaxy extend well beyond its disk. In a model with a 
fitted gaussian scale height h e ~ 20 kpc above the Galactic 
plane, and a scale radius p e ~ 35 kpc in directions perpen- 
dicular to the Galactic axis, we reproduced the observed 
properties of the GBR. To obtain these results, we as- 
sumed the nuclear CRs and the CR electrons to be equally 
distributed — as they would, if accelerated by CBs — , with 
their flux ratio fixed to the locally observed value. A scale 
height of CR nuclei as high as h CR = h e ~ 20 kpc rises the 
eyebrows of some experts and the swords of referees, but 
is not excluded by data on relative CR abundances 9 . A 
CR distribution of this large scale is a "cosmic-ray halo" . 
Its computed total luminosity, for the fitted values of h e 
and p e , coincides with the "large" value of Eq. QJ. 

A CR halo of similarly distributed electrons and ions is 
what is expected if CRs are accelerated by CBs, emitted 
by SNe close to the galactic disk. A CB of mass M CB 
and roughly constant cross section S CB , travelling in an 
interstellar medium of roughly constant density p ISM , has 
its Lorentz factor diminished from 70 to 7 in a distance: 



1 

7o 



(18) 



where = M CB /(p ISM S CB ). By explicitly fitting the 
spectra and time dependence of the afterglows of all GRBs 
of known redshift, Dado et al. (2002, 2003a) have found 
that the fitted values of distributed over two orders 

of magnitude, peaking at Xoo ~ 1 Mpc, while the values 
of 70 are narrowly peaked around 10 3 . It takes a distance 
£00/70 ~ 1 kpc for a CB to half its original Lorentz factor, 
while transferring half of its energy to CRs. We cannot 
extend with confidence Eq. I|18fl all the way to 7 ~ 1 
and trust the result: afterglow fluences decrease with a 
large power of 7 and cannot be followed long enough, and 
one does not expect the approximations of constant S CB 
and p 1SM to hold over distances much longer than a few 
kpc. Be it as it may, Eq. I|18fl indicates that it takes very 
much longer for a CB to lose and transfer the other half of 
its kinetic energy: 0(1) Mpc in the stated untrustworthy 
approximations. 



9 For the most elaborate models (Strong and Moskalenko 
1998) a "leaky-box" scale height of 20 kpc is only some 
1.3cr above the central value of the most precise observa- 
tions (Connell et al. 1998); it is perfectly compatible with 
the average of all previous and somewhat less precise results, 
l;2cl;2cl;2ccompiled in Lukasiak et al. 1994. 
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Table 1. In the list of clusters, Sersic stands for Sersic 159-03 
and EMSS for EMSS 1358+6425. Redshift 2, time-asymptotic 
central temperature Tf, inner observed temperature Tinner, 
outer observed temperature T ou ter, and input initial tempera- 
ture Ti, all in keV. The errors are given at the la level. 



45.0 
44.0 
43.0 
42.0 

46.0 
45.0 
44.0 
43.0 

42.0 
41.0 








b 


I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 

c 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ] 
d ! 



44.5 45.0 45.5 46.0 46.5 44.5 45.0 45.5 46.0 46.5 
log 10 L„ (ergs/s) log,,, L, (ergs/s) 



Fig. 1. The X-ray versus the optical luminosity of various sets 
of clusters, adopted from Miller and Nichol (2001). The thin 
(black) lines are best fits of the form Lx — a L" ; in the dashed 
lines a — 1 is enforced. The thick (red) line is the prediction of 
Eq. @, for which a ~ 0.02 and a = 1. (a) is the Abell/ACO 
sample, for which a = 0.90 ± 0.17. (b) is the EDCC sample, 
for which a = 0.95 ± 0.22. (c) is the CfA group sample, for 
which a = 1.06 ± 0.11. (d) is the RASS sample, for which 
a = 1.18 ± 0.08. 
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Fig. 2. a) A typical CF-cluster density profile, b) The time 
evolution of the temperature profile of the cluster, for the given 
density profile, a typical WR luminosity and an initial T = 3 
keV. The cluster ages correspond to 1/3, 1 and 3 times the 
current age of the Universe. The meaning of the four temper- 
atures, T, Tf, Ti n „er and T out er , is illustrated. 
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Fig. 3. a) The density profile of the cluster A1795. b) Its pre- Fig. 4. a) The density profile of the cluster Sersic 159-03. b) Its 
dieted temperature distribution. The data are from Tamura et predicted temperature distribution. The data are from Kaastra 
al. (2001). et al. (2001). The T outer value in Table I is that of the penul- 

timate radial point. 




Fig. 5. a) The density profile of the cluster Hydra, b) Its pre- Fig. 6. a) The density profile of the cluster A2029. b) Its pre- 
dicted temperature distribution. The data are from David et dieted temperature distribution. The data are from Lewis et 
al. (2001). al. (2002). 





Fig. 7. a) The density profile of the cluster A2052. b) Its pre- 
dicted temperature distribution. The data are from Blanton et 
al. (2002). The T ou t er value in Table I is that of the penultimate 
radial point. 



Fig. 8. a) The density profile of the cluster EMSS1358+6425. 
Its predicted temperature distribution. The data are from 
Arabadjis et al. (2002). 
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X-ray Mass [Solar mass] 

Fig. 9. a) The lensing and X-ray masses of a selection of clus- 
ters, from Hattori et al. (1999). The thinner (black) line is 
for equal masses and the wider (red) line is MLens = 3Mx, 
the result for "equipartition" of the plasma, the WRs and the 
magnetic field. 
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Fig. 10. Upper panel: The thick (red) line is a "x-by-eye" 
fit of Eq. 11141 to the thin thermal bremsstrahlung emission 
from the SN remnant SNR MSH 15-52. The data and figure 
are from Mineo et al. (2001). Lower panel: The same for the 
cluster A2256. The data and figure are from Fusco-Femiano et 
al. (2000). The thin (black) line is their binned, purely thermal 
extrapolation. 



